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Abstract 

We calculate supersymmetric two-loop corrections to the anomalous magnetic mo- 
ment of the muon, consisting of diagrams with a closed scalar fermion or fermion 
loop and gauge and/or Higgs boson exchange. We discuss the numerical impact of 
each subclass of diagrams and determine the leading contributions. We analyze in 
detail constraints from experimental information on the Higgs boson mass, Ap, and 
the branching ratios of i? — > Xgj and Bg ji'^ ji~ . If these constraints are taken into 
account, the largest possible effect of our two- loop corrections is reduced from more 
than 3 a (in terms of the current experimental error) to ~ 0.5 o", such that the influence 
on the total supersymmetric prediction is smaller than previously estimated. However, 
exceptions arise in rather extreme parameter scenarios with a strong non-universality 
between the soft breaking parameters in the stop and sbottom sectors. 
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1 Introduction 



A new era of precision measurements of the anomalous magnetic moment of the muon 
= {g^ — 2)/2 has been initiated by the "Muon g-2 Experiment" (E821) at BNL, leading 
to the current experimental world average of [1] 



a^^P = (11 659 208 ± 6) x 10"'" . (1) 

The most recent e^e~ data driven evaluations of the hadronic contributions by Refs. [2-4] 
lead to the following Standard Model (SM) predictions (the deviation from the experimental 
result is also shown) ^: 



^theo ^ (11 659 180.9 ± 8.0) X 10"^° (27.1 ± 10.0 
^theo ^ (11 659 175.6 ± 7.5) X 10-1° (32.4 ± 9.6 
^theo ^ (11 659 179.4 ± 9.3) X 10-1° (28.6 ±11.1 



2.7 a) [2] 
3.3 a) [3] 
2.5 a) [4] 



Recent analyses concerning r data indicate that uncertainties due to isospin breaking effects 
may have been underestimated earlier [4], so that with a better theoretical understanding of 
the isospin breaking effects the r-based results could come closer to the e+e-- based results. 
One may thus hope that eventually a combination of e"*'e~ and r data will lead to an even 
more precise theoretical prediction. 

While the present ~ 2.5 — 3.3(7 deviation between the SM prediction for and the 
experimental result can of course not be regarded as strong evidence for new physics, an 
increased accuracy of both theory and experiment might give rise to a significantly larger 
deviation in the future. On the other hand, already the current precision leads to very 
restrictive bounds on new physics scenarios. 

This is illustrated by the fact that the experimental precision of 6 x 10~^^ has now reached 
the level of the standard electroweak and of typical supersymmetric (SUSY) contributions. 
The electroweak one- and two-loop (and higher-order) contributions in the SM amount to 
19.5 X 10-1° and —4.1 x 10-i°, respectively, see Refs. [5,6] for reviews. The SUSY contribu- 
tions are generally suppressed by M^/M^, where Mw is the mass of the W boson and M is 
the typical scale of the SUSY particle masses. However, for large values of tan /?, the ratio 
of the vacuum expectation values of the two Higgs doublets in the Minimal Supersymmetric 
Standard Model (MSSM), the muon Yukawa couphng is enhanced by tan/3 as compared to 
the SM. The supersymmetric one-loop contribution is approximately given by [7] 



^ ' ^ M 



100 GeV\^ 



j tan/3, (2) 



where all SUSY masses are assumed to be equal to M. The involved SUSY particles are 
neutralinos, charginos and scalar leptons of the second generation. The magnitude of the 
supersymmetric one-loop contribution is at the right level to account for the ~ 3a deviation 
between the SM prediction and the data, and even larger shifts are possible. The supersym- 
metric two-loop contributions are known only in some approximations. Since the one-loop 

^The numbers for the combination of the experimental and the theory error and the corresponding 
deviaton in terms of a have been recalculated according to the new experimental result. 
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contribution can be large, the two-loop corrections can be expected to be quite important, 
even beyond the leading QED-logarithms [8]. 

The importance of the SUSY two-loop contributions is twofold. On the one hand, their 
inclusion increases the accuracy of the bounds on the supersymmetric parameter space (see 
e.g. Refs. [9, 10]). On the other hand, the supersymmetric two- loop contributions depend on 
many additional parameters and can in principle be large even if the one-loop diagrams are 
suppressed due to heavy smuons and sneutrinos. 

Particularly interesting contributions are the ones enhanced by large values of the Higgs 
mixing parameter fi and a large trilinear coupling A (where A generically denotes the Higgs- 
stop or Higgs-sbottom coupling, At^i,) . They arise from so-called Barr-Zee two-loop diagrams 
where a Higgs boson is exchanged between the external muon and a 3rd generation sfermion 
loop. Results for such contributions were obtained in Refs. [11,12] and found to give huge 
contributions up to 0(20 x 10~^°) if tan f] is large and /i, A are of the order of several TeV. In 
these analyses, however, other experimental constraints on the parameter space of the MSSM 
were neglected. Moreover, the results of Refs. [11, 12] for the contribution disagree by a 
factor 4, so that an independent check seems to be necessary. 

In this paper we present a calculation and numerical analysis of all two-loop contributions 
Aa^'" in the MSSM where a closed 3rd generation sfermion or fermion loop is inserted into 
a one-loop diagram with gauge-boson and/or Higgs-boson exchange. This set of diagrams 
contains the terms oc /i, A but also other terms enhanced by the large Yukawa couplings of 
the t, and (for large tan/3) b, r, as well as terms without any enhancement. All of these 
contributions are included in our final result. 

Our numerical analysis is focused on two questions: what are the numerical results 
for the individual subclasses, and which of them should be taken into account for a reliable 
supersymmetric prediction for a^? Secondly, are huge two-loop contributions of C(20x 10~^°) 
still possible if existing experimental constraints on the supersymmetric parameter space are 
taken into account? 

The rest of the paper is organized as follows. In Sect. El we present the two-loop dia- 
grams and the method of their evaluation. The numerical analysis is given in Sects. El S El 
The importance of this class of two-loop corrections and their numerical size is discussed in 
Sect. El taking into account constraints on the SUSY parameter space from other experimen- 
tal information. The leading contributions and the influence of the individual experimental 
constraints are examined in Sect. |31 Effects from non- universality of the soft SUSY-breaking 
parameters are analyzed in Sect. El We conclude with Sect. El 

2 Calculation 

In this section we briefly describe the diagrams we have investigated, their evaluation and 
the tools that have been used. 

The set of diagrams calculated in this paper corresponds to the fermion/sfermion correc- 
tions to non-supersymmetric, i.e. two Higgs doublet model type, contributions. It forms a 
gauge-independent class of diagrams. In order to discuss the shift between the MSSM and 
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the SM predictions, we subtract the pure SM contribution from our result (where the SM 
Higgs boson mass Mf^ is set to the value of the lightest MSSM Higgs boson mass, M/j). 

The one-loop diagrams corresponding to the contributions considered in this paper are 
the SM one-loop diagrams. Expressing the one-loop result in terms of the Fermi constant 
Gfj, and = 1 — My^/M'^, {Mz being the Z boson mass), it takes the conventional form of 
the electroweak one- loop result in the SM (omitting the QED contribution) [5,6], 



RW,1L ^ _^^2 



5 + i(l-4.?.r 



(3) 



The two-loop diagrams that we calculate can be subdivided into three classes: 
ifVcp) diagrams with a sfermion (t, b, f, z/^) loop, where at least one gauge and one Higgs 
boson are exchanged, see Fig. HJ 

(fVV) diagrams with a sfermion loop, where only gauge bosons appear in the second loop, 
see Fig. |21 

(fVcp) diagrams with a fermion {t, b, r, z/,-) loop, where at least one gauge and one Higgs 
boson are present in the other loop, see Fig. El The corresponding diagrams with only gauge 
bosons are identical to the SM diagrams and give no genuine SUSY contribution. 

For our later analysis we further split up the {fV(f)) diagrams of Fig. ^ into the fol- 
lowing groups: diagrams with photon and Higgs exchange {f'y{h,H}), Z/Higgs exchange 
{fZ{h,H}), and W^/charged Higgs exchange {fW^H"^). The remaining sfermion loop di- 
agrams containing only gauge boson and Goldstone boson exchange are grouped together 
with the diagrams involving only gauge bosons, {fW'^G^) + {fVV). All these groups are 
separately gauge independent in _Rg-gauges at the order m^/M^. Note that since we neglect 
CP- violating phases, diagrams with photon or Z and CP-odd Higgs bosons A'^, do not 
contribute. The diagrams {f'j{h, H}) and (fW^H^) are the ones evaluated in Refs. [11, 12] 
neglecting all but the leading terms in the sfermion-Higgs couplings. 

All diagrams are understood to include the corresponding subloop renormalization. For 
the fermion loop class (fVcj)) we actually calculate the difference between the Standard Model 
and the MSSM, which originates from the extended Higgs sector of the MSSM. Diagrams 
where two Higgs bosons couple to the external muon are suppressed by an extra factor of 
m^/M^ and hence negligible. 

In order to perform a systematic calculation, a Mathematica program has been written 
that can deal with all kinds of MSSM two-loop contributions to a^. Its main steps are the 
following: The amplitudes for are generated using the program FeynArts [13,14], and the 
appropriate projector [15, 16] is applied. The Dirac algebra and the conversion to a linear 
combination of two-loop integrals is performed using TwoCalc [17]. In order to simplify the 
integrals, a large mass expansion [18] is applied where the muon mass is taken as small and all 
other masses as large. All resulting two-loop integrals are either two-loop vacuum integrals 
or products of one-loop integrals. They can be reduced to the standard integrals T134 [19] 
and Aq and Bq [20] and can be evaluated analytically. The asymptotic expansion has to 
be performed up to terms of order m^. Terms of lower power in m^/Mheavy (where Mheavy 
represents all kinds of other masses) cancel each other as required for non-QED corrections. 
Terms of higher powers in ?7i^/Mheavy are numerically irrelevant and can be safely neglected. 

The counterterm diagrams contain the renormalization constants 6M^^, 6Ze, 6th^H cor- 
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Figure 1: Some generic two- loop SUSY diagrams of type {fVcp) involving (at least) one gauge 
and one Higgs boson and a closed scalar fermion loop. F = /x, z/^; (f) = h,H,A, H^, G, G^; 
■0 = G^; f,f' — t,h,f,v-,-] V — ^,Z,W. Further diagrams of this type are obtained by 
contracting the line in the first diagram or by interchanging Higgs and vector bosons. 




f 
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Figure 2: Some generic two- loop SUSY diagrams of type (fVV) involving gauge bosons and 
a closed scalar fermion loop. F — /i, u^; f, f = t, b, f, z/^-; V = •y, Z, W. Further diagrams of 
this typ involving four-point vertices exist as well. 




Figure 3: Generic two-loop SUSY diagrams of type {fV(j)) involving (at least) one gauge and 
one Higgs boson and a closed SM fermion loop. F = /i, u^; 4> = h, H, A, H"^, G, G^; ip — G^; 

f,f' = t,b,T,Ur;V = ^,Z,W. 
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responding to mass, charge and tadpole renormalization and can be easily evaluated. We 
choose the on-shell renormalization scheme [21]. This leads to 5M^^^ = ReE^^2(M^^^), 
where S^^^ denote the transverse parts of the gauge-boson self-energies. The charge renor- 
malization is given by ^Z^ = —1/2 T,'^{0), where S' denotes the derivative of the self-energy 
with respect to the momentum squared. The tadpoles are renormalized such that the sum 
of the tadpole contribution T and the counterterm vanishes, i.e. Sth^n = —Th,H- 

As mentioned above, see eq. 0, we are using a one- loop result which is parametrized 
in terms of instead of the ratio a/M^. Therefore our two-loop correction contains a 
term given by the product of the corresponding one-loop result and Ar, where the latter 
denote the one-loop corrections to muon decay, fi — *• v^eve- Relevant here are only the 
contributions arising from 3rd family sfermion loops to Ar. These corrections are included 
in our two-loop result (in the {fVV) class). 

As a cross check we have evaluated the SM two-loop diagrams with a closed fermion 
loop as presented in Refs. [15, 16] and found perfect agreement separately for each diagram 
(after going to the limit — 1/4, used in Refs. [15, 16]). We have furthermore checked the 
UV-finiteness of our result as well as the cancellation of the wave function renormalization 
constants. We also found agreement with Ref. [12] for the contributions to the {f'y{h, H}) 
and {fW^H"^) diagrams calculated there. This confirms that the earlier result of Ref. [11] 
is too large by a factor of 4. 

Our final result for the sum of all diagrams is rather lengthy and not displayed here. It 
is included as a Fortran subroutine in the code FeynHiggs [22] (see: www.feynhiggs.de). It 
can also be obtained as a Mathematica formula from the authors upon request. 



3 Numerical results allowed by experimental 
constraints for the different sets of diagrams 

The MSSM two-loop contributions to depend on many parameters, most notably on tan (3, 
the fi parameter, the trilinear soft SUSY-breaking parameters At^^r, the mass of the CP-odd 
Higgs Ma, and the soft SUSY-breaking parameters Mq^l^u,d,e appearing in the sfermion 
mass matrices. In Refs. [11,12] it was shown that in particular large /i and A parameters 
can give rise to very large contributions of the {/^{H, h}) and [fW^H"^) diagrams, however 
ignoring existing experimental constraints on the MSSM parameter space. In order to find 
out the largest possible contributions of each class of diagrams, we perform a scan of the 
MSSM parameter space. We vary the parameters in the ranges 



-3 TeV 


< 


/i 


< 


3 TeV 


-3 TeV 


< 


At,b 


< 


3 TeV 


150 GeV 


< 


Ma 


< 


1 TeV 





< 


^SUSY 


< 


1 TeV 


Mq = Ml 




Mu = 


Mn 


= Me 



(4) 

where we have set Msusy 
breaking parameters in the t mass matrix, Mq,Mu in the b mass matrix, and Ml,Me in 
the f mass matrix) and Aj. = Ai,. Furthermore we fix tan/? to tan/3 = 50. Large values 
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Quantity 


Mh 




BR(5, ^ 11+ fi-) 




strong bound 
weak bound 


> 111.4 GeV 

> 106.4 GeV 


< 3 X 10^^" 

< 4 X 10-1° 


< 0.97 X 10-^ 

< 1.2 X 10-^ 


< 1.0 X 10"^ 

< 1.5 X 10"^ 



Table 1: Strong and weak bounds imposed on the MSSM parameter space. As^Xs-y = 
|BR(i? — > 3.34 X 10-^1, where 3.34x 10"^ is the current experimental central value [32]. 

of tan/? and small values of M4 generically lead to larger SUSY contributions but also 
to more restrictive experimental constraints. These two effects tend to cancel each other. 
We have checked that the maximum contributions from our diagrams to allowed by the 
experimental constraints are about the same for tan j3 = 25, tan f3 = 37 and tan (3 = 50 and 
when Ma is varied in the range Ma = 90 . . . 150 GeV. The effect of relaxing the restriction 
of a common soft SUSY-breaking parameter in the sfermion mass matrices will be described 
in Sect. El As SM input parameters we use rrit = 175 GeV, and mb{mt) = 3 GeV (in order 
to absorb leading QCD corrections). 

We restrict the parameter space further by imposing the following experimental con- 
straints:^ 

• The lightest CP-even MSSM Higgs-boson mass Mh has to be larger than its experi- 
mental limit 114.4 GeV [25,26].'^ has been evaluated with FeynHiggs2.0 [22], based 
on Refs. [27-29]. 

• The t/6-contribution to the p parameter, evaluated up to the two- loop level [30], does 
not exceed its experimental bound. 

• The branching ratios BR(i?s P^f^^) [31] and BR(i? Xg'j) [32] are in agreement 
with their experimental limits.^ 

In order to be able to check the sensitivity on these bounds we use a stronger and a weaker 
version for each bound, see Tab. [T] The two Higgs mass bounds take into account a 3 GeV 
uncertainty due to unknown higher-order corrections [29], the weak bound in addition an 
uncertainty of 5 GeV due to the imperfect knowledge of the top mass [33], on which is 
much more sensitive than a^. The two bounds on Ap^^^^, BR(i? —>■ Xg'-f), and BR(i?s — > 
fi+fi~) correspond to 2a and 3a bounds and to 90% and 95% C.L. bounds, respectively. 

The interplay of these constraints restricts the allowed parameter space severely. The 
Mh bound puts a limit of about 2.5Msusy on \At\. The data on b decays constrain /i and 
A parameters in particular for small Ma- Ap restricts the mass splittings in the t, b sectors 

^There are of course also lower bounds on sfermion masses from direct searches. The bounds from LEP 
are roughly rrif i ^ 100 GeV. From Run I of the Tevatron stronger bounds arise for parts of the MSSM 
parameter space [23,24]. We do not impose the direct bounds explicitly in the scans since we present the 
results of as functions of the lightest sfermion mass. 

^The limit on the SM Higgs mass holds unchanged for the light MSSM Higgs mass Mh for Ma Z 150 GeV. 
For lower values of Ma , the bound on Mh is smaller but the restrictions implied on the and A parameters 
are significant also in this case. 

^We are grateful to A. Dedes and G. Hiller for providing the respective codes. 
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Figure 4: Possible contributions to Aa^^ for the case that all experimental bounds are 
required in their strong versions. The results are subdivided into five classes of diagrams: 
sfermion loops with gauge and Higgs boson exchange {f'y{h,H}), (fW^H^), {fZ{h,H}), 
sfermion loops with gauge or Goldstone boson exchange {fVV) + {fW'^G^), and fermion 
loop diagrams {fV(f)). The results are plotted as functions of the lightest sfermion mass 
(sfermion loops) or Ma (fermion loops). 
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and thereby also the /i and A parameters, which appear in the off-diagonal elements of the 
squark mass matrices. 

In Fig. m we plot the resulting ranges of possible contributions of the individual classes 
of diagrams for the case that all bounds are required in their strong versions. In the case 
of the sfermion loop contributions we plot the results for over the lightest sfermion mass 
{mm{mf^,mt^,mi_^,mi^}), and in the case of the fermion loop contributions we plot the 
results over Ma- 

From Fig. |3] the following conclusions can be drawn: 

• The contribution of the gauge and charged Goldstone boson exchange diagrams, 
(fW^G^) + ifVV), is very small. Its maximum size is about ^af^'^^ ^ 0.02 x 10"^°. 

• The contribution of the {fZ{h, H}) diagrams with Z and Higgs exchange is at most of 
the order 0.1 x 10^^° and thus negligible compared to the present experimental error. 
The reason for this suppression compared to the photon and W exchange diagrams is 
the factor (1 — 4s^) in the coupling of the Z to muons. 

• The contribution of the {fV(f)) diagrams with a fermion loop can reach 0.6 x 10~^° for 
small Ma ^ 200 GeV; the fermion loop diagrams are thus not completely negligible. 

• The photon exchange diagrams with a sfermion loop (/70) are dominant; the results of 
the W exchange diagrams (fW^H^) are much smaller. The photon exchange diagrams 
are the only ones that can contribute more than 1 x 10~^°, the W exchange diagrams 
contribute up to 0.3 x 10"^''. The reason for the suppression of the W diagrams is 
not only the high value of M\y but also the fact that the W couples to two different 
sfermions, to t — 6 or to f — !>,-, of which at least one is usually relatively heavy. 

Hence we find that the photon exchange contributions calculated in Refs. [12] are indeed 
the dominant subclass of diagrams with a closed (s)fermion loop. We also find, however, 
that the maximum contributions of more than 20 x 10~^° quoted in Refs. [11,12] for the 
photon and W contributions are reduced to about 2.5 x 10~^° and 0.3 x 10~^° due to the 
experimental constraints on the MSSM parameter space. Owing to the smallness of these 
contributions, the fermion loop contributions can make up for a non-negligible part of the 
two-loop corrections. 

4 Leading contributions and influence of the experi- 
mental constraints 

Let us now focus on the photon and W exchange contributions {fj{h,H}), {fW^H"^) and 
study the influence of the individual constraints on the ranges of possible numerical values. 
We choose this set of contributions not only because the photon contributions are dominant 
within our class of diagrams, but also because these contributions are obviously significantly 
restricted by the experimental constraints. In contrast, the fermion loop contributions (/V^0) 
can be non-negligible but depend mainly on Ma and are hardly constrained. 
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It is instructive to explicitly discuss the complete expression for the photon diagrams 
(see e.g. Ref. [12]): 



2 



where / can be one of ii^2, ^1,2, ^1,2, and (f) can be one of the CT'-even Higgs bosons, h or H. 
The couplings A are defined as = sin a , Ca — cos a , etc.) 

^^i{h,H} = {-Sa,Ca}/Cf3 (6) 

>^ii{h,H} = 2mt(/i{s„,-Ca} + At{c„,s«})f7^iL^2/s/3 



(t^^)' (7) 



w /TTt ^2 



= 2mb( -/i{c„,s«} + A{-s«,c„})f/fif/f2/c/3 

^ 6cvymg{-Sa, Ca} + MwMzCf}{-?, + 2g^){-ga+;3 , Ca+j3 } ^^^^ 2 



iU^2? (8) 



%Cwml{-Sa, Ca] " 2{-Sa+l} , C^+p }MwMzCps'^ ^^^fo ^2 



and similar for \f{h,H}- The matrices ?7*'^ diagonalize the sfermion mass matrices M| in the 
form M'jiU^y = diag(m^^, m|^). The loop function T is given by 

^ ri^^ a;(l-x)log[./(x(l-x))] _ 
Jo 2; — x(l — x) 

The result for the W contribution has a similar form. 

This type of contributions can be particularly enhanced by the ratio of the mass scale of 

the dimensionful Higgs-Sfermion coupling divided by the mass scale of the particles running 

2 

in the loop, i.e. by ratios of the form {/i, >1, -^j}/{mj, M/j^i^}, which can be much larger 
than one. For large tan /3 and large sfermion mixing, the leading terms are typically given by 
the parts of the couplings with the highest power of tan f3 and by the loop with the lightest 
sfermion. These contributions involve only iJ-exchange, since the /i-couplings approach the 
SM-Higgs coupling for not too small M^. They can be very well approximated by the 
formulas 

Aa*/L = _o.013 X 10-1° ^^^^^^^ sign(A), (10) 

miMH 

Aa^'^L = -0.0032 X 10-1° ^^^^^^^ sign(/.), (11) 

where mi and are the masses of the lighter i and 6, respectively, and Mh is the mass of 
the heavy CP-even Higgs boson. The formulas use the approximation J^{m'j/M'^)/M]j ft; 



0.34/ {mjMH) for the loop function, which holds up to few percent if the respective sfermion 
mass fulfils mn ^ M^. Since the heavier sfermions also contribute and tend to cancel the 
contributions of the lighter sfermions, these formulas do not approximate the full result very 
precisely, but they do provide the right sign and order of magnitude. 

Equations (fTUj) . show that the mj-contributions are enhanced by one power of tan/5 
from the muon Yukawa coupling and by the ratio fi/Mn, whereas the mt-contributions 
contain an additional power of tan/3 from the b Yukawa coupling and the ratio Ab/Mn- For 
tan/3 = 50 and /z and A parameters larger than 1 TeV both contributions can amount to 
more than 1 x 10~^°. However, /i is much more constrained by the four experimental bounds 
in Tab. H than A^. Therefore, the largest contributions in Fig. |3] originate from the sbottom 
loop diagrams and from parameter constellations where bi is the lightest sfermion. 
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Figure 5: Maximum contributions of the (/7{/i, -f^}) and (/VT if^) diagrams to Aa^ as a 
function of the lightest squark mass, min{mj^, m^^, "m^^, mi^. No constraints except for the 
parameter ranges in eq. (@]) are taken into account for the outermost curve. Going to the 
inner curves additional weak constraints (see text) have been applied. 



Now we study the influence of the individual experimental constraints on the photon 
and W exchange contributions. Figure El is based on a data sample of ~ 300000 parameter 
points in the range specified in eq. (j^, on which the weak versions of the bounds in Tab. 
are incrementally applied. FigurelHlis based on the data points satisfying all weak constraints 
and shows the effect of strengthening each bound separately. 

The results shown in Fig. are the following: 

• The outer lines show the largest possible results if all experimental bounds are ignored. 
They show a steep rise of Aa^'" for decreasing m^^; for < 150 GeV contributions 
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50 100 150 200 250 300 350 400 450 500 
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Figure 6: Maximum contributions of the {f^{h, H}) and (fW^H^) diagrams to Aa^^ as a 
function of the hghtest squark mass, minjm^^, mf^, mi_^, 'm-i^}. The outer curve corresponds 
to weak bounds for all experimental constraints. Each inner curve takes into account one 
additional stronger constraint. Strengthening the Bs —>■ /x"'"/x~-bound has a very small impact 
and is not shown. The inner area consequently corresponds to all strong constraints. 



larger than 15 x 10 corresponding to two standard deviations of the experimental 
error on a^, are possible. 

• The next two lines show the possible results if the bound > 106.4 GeV and then 
in addition the bound on Ap are satisfied. The maximum contributions are very much 
reduced already by the Mh bound, and the Ap bound reduces further the positive region 
for small sfermion masses. If both bounds are taken into account, Aa^^ > 5 x 10~^° 
and Aa^L < -10 x 10-^° is excluded for mj^ > 100 GeV. 

• The two innermost lines correspond to taking into account in addition the bound on 
BR(i?s H.'^IjT) and finally also on BR(S — > Xg'^). In particular taking into account 
the BR(S — > Xg^) bound ehminates most data points with mj^ < 150 GeV and 
thus leads to a strong reduction of the possible size of the contributions. The largest 
contributions of ±4 x 10^^*^ to Aa^'", corresponding to ~ 0.7cr of the experimental 
error, are possible for rrij^K, 150 . . . 200 GeV. 

In applying all bounds one should be aware that any fiavour non-universality in the MSSM 
parameters could have a strong effect on the predictions for the h decays, whereas the influ- 
ence on Mh and Ap^^^^ would be mild. Hence it is interesting that even if the h physics 
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bounds are ignored and only the weak and Ap bounds are taken into account, the largest 
possible contributions are strongly restricted to 

- 10 X 10"^° < /\af < 5 X 10"^° (12) 

for sfermions heavier than 100 GeV. 

Figure IHl shows that strengthening the bound on from > 106.4 GeV to > 
111.4 GeV has the most significant effect. It cuts off all the regions where Aa^'" > 3 x lO"^'' 
and Aa^^ < —2 x 10~^°. Strengthening the other bounds has only a marginal effect. This 
confirms that the M/j-bound is most important for restricting the parameter space. 



5 Non-universal soft SUSY-breaking parameters 

Up to now we have found only moderate numerical effects from the two-loop diagrams with a 
closed (s)fermion loop, even for the photon exchange diagrams. However, the approximation 
formula eq. (fTUj) . Aa^^'" oc fimt/^m^MH), shows that values up to 15 x 10^^° should be 
possible if /i ~ 3 TeV and m^-, Mh ~ 150 GeV (for tan/? = 50). In Fig. El such contributions 
indeed appear but they are excluded if the experimental constraints are taken into account. 
Already imposing only the bound on reduces the maximum contributions almost by a 
factor of 3. 

It is important to keep in mind that all results presented so far were based on the 
universality assumption Mq = Mu = Mr, = Ml = Me = Msusy for the soft SUSY-breaking 
parameters in the MSSM.^ In this section we examine the effect of relaxing this assumption. 
We do not attempt a full scan of the MSSM parameter space but rather investigate which 
pattern of non-universality can lead to particularly large results. 

The reason why the large results in Fig. are excluded is that universality indirectly 
leads to severe constraints on the fi parameter. Via universality, the left- and right-handed 
diagonal elements of the stop and sbottom mass matrices are linked, and if one requires 
a light stop, there is not much room for an even lighter sbottom. Hence the off-diagonal 
element in the sbottom sector cannot be much larger than the one in the stop sector, which 
means for large tan/?, tan (3 > rrit/mi,: 

l/^l ^ l^tl- (13) 

At is not only restricted by the requirement that all stop squared masses are positive but also 
by the M/i-bound, which roughly leads to \At\ < 2.5Msusy- Light stops require [mfAtl ^ 
MgusY and are therefore only possible for Msusy ^ 400 GeV, and thus \At\ can hardly 
exceed 1 TeV. Because of universality and eq. (fT^ . the bounds on At hold also for /x, and 
therefore also /i ^ 1 TeV. 

^SU(2) gauge invariance dictates only that the same left-handed squark mass parameter Mq appears 
in the stop and sbottom mass matrices (and analogously for Ml in the slepton sector). Apart ft'om this, 
there is a priori no reason in the unconstrained MSSM to assume equality of the left- and right-handed 
sfermion mass parameters, except for simplicity. The symmetries of the MSSM allow independent values for 
Mq,u^d,l,e- 
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8000 



= 400 GeV 





Figure 7: Possible values of fi (upper plot) and corresponding contributions to Aa^'" (lower 
plot) for the case of non-universality of the soft SUSY-breaking parameters. The plots show 
^ and Aa^^ as a function of Msusy = Mq = Mjj for different values of the ratio M/^/Msusy- 
The other parameters are chosen as m^^ = 150 GeV, M4 = 400 GeV, A^, = 0, tan/5 = 50. 



It is therefore interesting to break up the relation between the stop and sbottom mass 
parameters and to require only 

Msusy = Mq = Mu = Ml ^ Md = Me. (14) 

Thus we can choose small values of Msusy and At, giving rise to a light stop. Choosing 
Md S> Msusy at the same time allows very large fj, without producing a too light sbottom. 
We will see that these large values of fj, are also compatible with the bound on M^. 

Fig. [7| shows the results for different ratios of M/j /Msusy- We choose a light stop mass 
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r/ifj = 150 GeV and a moderate value Ma = 400 GeV in order to avoid too strong restrictions 
from h decays. For each Msusy, is determined by m^^ = 150 GeV. The values of /i are 
determined as the maximum values compatible with > 111.4 GeV and Ap^^^^ < 0.004.^ 

The upper plot in Fig. [3 shows these maximum values of /i as functions of Msusy- They 
significantly increase with M^/Msusy- Already for M^/Msusy = 3, values for larger than 
1.5 TeV are possible. For Md/Msvsy = 6, = 3 TeV is possible, and for Md/Msvsy = 30, 
even fi = 6 TeV is possible. 

The lower plot in Fig. [3 shows the corresponding results of the photon exchange diagrams 
Aa^J"'^'^'^^^''^^. We choose At < so that the contribution to is positive. The results 
exhibit a clear correlation with the values of fi, and they are quite precisely given by the 
approximation pO|) .'' Thus the maximum results with /i < 3 TeV are about 10 x 10~^°, and 
the results for /i = 6 TeV are larger than 20 x 10~^°. 

It should be noted that these parameter choices are rather extreme and involve vastly 
different mass scales for the MSSM parameters. As an example, the largest results are 
obtained for Msusy ~ 300 GeV, ~ 9 TeV, At ~ 550 GeV, /i ~ 6 TeV. We have checked 
that all points plotted in Fig. [3 with Msusy = 260 . . . 390 GeV, where fi and Aa^ are large, 
satisfy not only the bounds on Mh and Ap but also those on _B — Xgj and Bg 
Only if smaller values. Ma < 400 GeV, are chosen, strong violations of the b decay bounds 
occur. For larger Ma, on the other hand, the b decay constraints are less restrictive, and 
even larger values for p, and AaJ}" than in Fig. [7| are possible. 



6 Conclusions 

We have obtained results for MSSM two-loop corrections to the anomalous magnetic mo- 
ment of the muon. The corrections consist of diagrams where a SM fermion or sfermion 
loop is inserted into a one-loop diagram with gauge- and/or Higgs-boson exchange. We 
have investigated the importance of the individual contributions and the impact of existing 
experimental constraints on the maximum numerical results. 

It has been found that the by far most important of the considered diagrams are the ones 
with a sfermion loop and photon and neutral Higgs exchange {/^{h, H}). They contribute 
up to about 2.5 x 10"^" in the parameter space allowed by all experimental constraints. This 
value has to be compared with the current experimental error of 6 x 10^^*^. The diagrams 
with sfermion loop and W^/H"^ exchange (fW^H^) and the fermion loop diagrams {fV(f)) 
contribute up to 0.3 x 10~^° and 0.6 x 10~^°, respectively, while the remaining diagrams are 
negligible. 

Our second result is that taking into account existing experimental constraints is crucial. 
We have carefully analyzed the impact of the constraints on the lightest Higgs-boson mass, 
Ap, BR(i?s — >■ and BR.{B Xg'y). If the experimental constraints were ignored 

and the p, and A parameters were varied up to 3 TeV, contributions of more than 15 x 

^Ab is set to zero here since the sbottom contributions cannot be expected to increase significantly beyond 
~ 5 X 10-1°, see eq. (tTTll . 

^Owing to the large values of /i the loop corrections to the heavy CP-even Higgs mass Mh can be large, 
and Mh, which enters in eq. (|10|) . can be significantly lower than Ma- 
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10"^", corresponding to 2.5a of the experimental error, would be possible from the two-loop 
diagrams. Already if only the experimental bounds on and Ap are taken into account, 
the accessible parameter space for fi, At, and Aj, is severely restricted, and one obtains 
—10 X 10~^° < Aa^'" < 5 X 10~^°. Taking into account all constraints leads to the relatively 
small result of Aa^'" ^ 3 x 10^^°. This two-loop correction of ~ 0.5 o" therefore gives rise to 
only a moderate shift of the one-loop SUSY result (which can easily account for the ~ 3 a 
deviation between the SM prediction and the data). 

The results quoted above have been obtained under the assumption of universal soft 
SUSY-breaking parameters Mq = Mu = Mj:,. If one allows large mass splittings between 
these parameters, the considered MSSM two-loop contributions can have a significantly 
larger numerical effect while the existing constraints are still satisfied. We have analyzed 
the example of Mo > Mq = Mjj, which can give rise to particularly large contributions to 
a^. One needs large ratios Md/Mq > 3 and at the same time a light stop and extremely 
large fi in order to obtain contributions that are significantly higher than 5 x 10~^°. Though 
in principle possible, such parameter constellations look quite artificial, and they should be 
viewed as an illustration of how difficult it is to produce larger contributions to a^. Models 
with universality at some high scale typically lead to approximate low-energy universality 
Mq ^ Mu ^ Mu and Mq ^ M^ ^ Me, which would restrict the allowed range for fi even 
more than low-energy universality. 

The contributions presented in this paper involve the potentially large enhancement 
factors {/i, A}/{mj, Mh} and constitute an important part of the two-loop contributions in 
the MSSM. Our full result is included as a Fortran subroutine in the code FeynHiggs (see: 
www.feynhiggs.de). It can also be obtained as a Mathematica formula from the authors 
upon request. 

In order to reduce the theoretical uncertainty of the MSSM prediction for further, 
the remaining two-loop contributions should be analyzed as well. The technical tools devel- 
oped in this paper allow such a study, and the results will be presented in a forthcoming 
publication. 
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